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I ■ Abstract. The phenomenology of LiFeAs superconductor can be explained in the 

O ,, framework of a four-band s±-wave Eliashberg theory. We have examined the 

experimental data available in literature and we have found out that it is possible 
to reproduce the experimental critical temperature, the gap values and the upper 
critical magnetic field within an effective model in moderate strong-coupling regime 
that must include both an intraband term An ~ 0.9 and an interband spin-fluctuations 
1 . , (Kit ~ 1-5) coupling. The presence of a non negligible intraband coupling can be a 

^ | fictitious effect of the violation of Migdal's Theorem. 
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The family of iron pnictide superconductors, discovered by the Hosono group [T] 
in 2008, has focused an intense research in this four years. Many compounds with 
different crystal structures and physical properties have been discovered and character- 
ized [21 |3l H], among them LiFeAs [5] has proved to be a peculiar one. 
First of all LiFeAs, unlike almost all other iron- compounds, does not need neither charge 
doping nor pressure to condense in the superconducting state |5], this implies that no 
disorder is present. Further, it seems to be not magnetic and angle-resolved photoemis- 
sion spectroscopy (ARPES) reports poor nesting |6J. At a glance these characteristics 
could turn us away from the idea of an unconventional pairing mechanism, however the 
phonon-mediated coupling seems to be too weak [7] to explain the relatively high criti- 
cal temperature (T c =18 K). Moreover, recent neutron inelastic scattering measurements 
show the evidence of strong antiferromagnetic fluctuations [8], reconciling LiFeAs with 
other Fe-based superconductors. 

A multigap scenario is suggested by the presence of five different sheets in the Fermi 
surface [9]: Two electron pockets are centered near the M-point of the Brillouin zone 
and three hole pockets around the T-point. Despite the Fermi surface shows five 
different sheets, according to our electronic structure calculations the 5-th sheet can 
be disregarded because of its low density of states (see Table [1] and Figure [TJ and 
size pm [HI [121 1131 ITS] - Consequently, as a starting point, we can model the electronic 
structure of LiFeAs by using a four-band model [151 02] with two hole bands (1 and 2) 
and two electron bands (3 and 4). 

In this work we construct a theoretical model to describe the phenomenology of LiFeAs, 
by using a minimal number of phenomenological parameters in combination with Density 
Functional Theory (DFT) calculations and the Eliashberg theory of superconductivity. 
The model is then tested by simulating the temperature dependance of the critical field. 

Recently, four different gaps have been observed by ARPES measurements on 
LiFeAs [15] . Hence, in order to describe the superconductive properties of this 
compound, a four-band Eliashberg model [151 02] with s± symmetry [17J can be used. 
The experimental gap values [18] have been used to fix the free parameters of the model 
and, after this, the critical temperature and the upper critical magnetic field [19] have 
been calculated. The final result is a moderate strong-coupling regime Xtot ~ 1-6 — 2.0 
where the total electron-boson coupling must include two different contributions: a 
purely interband coupling mediated by spin-fluctuations (s/), with A^ t ~ 1.5 and a 
purely intraband coupling An, whose origin will be discussed hereafter. The isotropic 
values of the gaps at T = 8 K are reported to be Ai = 5.0 meV, A 2 = 2.6 meV, 
A 3 = 3.6 meV, A 4 = 2.9 meV. As a first approximation, since just a small anisotropy 
is observed, we consider only the isotropic part. 

The Eliashberg theory [20] generalized to multi-band systems has already been 
successfully used to describe the properties of MgB 2 [211 1221 [23] and iron compounds 
[24"1 [25] . A four-band Eliashberg model includes eight coupled equations for the gaps 
Ai(iu n ) and the renormalization functions Zi(iu n ), where i is the band index (that 
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Figure 1. Fermi Surface of LiFeAs. 



Table 1. Fermi Surface resolved Kohn Sham properties 10 : The Fermi density of 
states (iV(0)) is given in states/spin/eV, the Fermi Velocities (vf) in 10 5 m/sec, and 
plasma frequencies (w p ) in eV. ab label the in-plane and c for the out-of-plane direction 
of the Fermi velocities and the diagonals of the plasma tensor [14] . 
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ranges between 1 and 4) and u n are the Matsubara frequencies. The imaginary- axis 
equations are: 



u n Zi(iu n ) =u n + 7iT y^ AfjjiWn, ioj m )Nf{iu m ) + 

J 

Zi(iu n )Ai(iu n ) = nT^2[A^(iu n ,iu m ) - A*y(w c )] X 

x q(u c - \u m \)N*{iu m ) + £[r y + rf ]jv/ 



(i) 



(lUr, 



(2) 



where T;j and ly are the non magnetic and magnetic impurity scattering rates, 
Afjiiu^iUm) = A p ^(iu n ,iu m ) + A$ {iu n , iu m ) , A§(iu n ,iu m ) = A p ^(iu n ,iu m ) - 
A%jj(iu) n ,iw m ), sf means antiferromagnetic spin fluctuations and ph phonons. 
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6(cj c — \u> m \) is the Heaviside function and u c is a cutoff energy. In particular, 
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Coulomb pseudopotential matrix and, finally, 
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The electron-boson coupling constants are defined as 
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nevertheless. 



a huge number of 
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The solution of Equation ([I]) and Equation 
input parameters (32 functions and 16 constants) 
interdependent, others may be extracted from experiments and still others fixed by 
appropriate approximations. 

At the beginning we fixed the same conditions that have been used for many other 
pnictides, as shown in Ref. [17] . and we assumed that: i) the total electron-phonon 
coupling constant is small and mostly intraband [26]; ii) antiferromagnetic spin 
fluctuations mainly provide interband coupling [31 [24] . To account for these assumptions 
in the simplest way (as has already been done for other iron- compounds with good 
results) we should take: A 



ph 
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ph 

ij 



0, /i*j(w c ) = lAjd^c) = i.e. the electron-phonon 
coupling constant and the Coulomb pseudopotential compensate each other, in first 
approximation, and X S J = 0, i.e. SF produce only interband coupling [24J. However, 
within these assumptions, we were not able to reproduce the gap values of LiFeAs, and 
in particular the high value of Ai, the best results obtained are reported in the second 
row of Table |2j 

In order to solve this problem it is necessary to introduce at least an intraband coupling 
in the first band, then An ^ 0. 
The final matrix of the electron-boson coupling constants becomes 
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where z/jj = Ni (0)/Nj(0) and Aj(0) is the normal density of states at the Fermi level for 
the i — th band (i = 1, 2, 3, 4). We chose spectral functions with Lorentzian shape i.e: 
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and Cij are normalization constants, necessary 
to obtain the proper values of \j while Qij and Y^ are the peak energies and half- 
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because the ARPES measurement are on very good single crystals of LiFeAs [18J. Band- 
structure calculations (see Table [p provide information about the factors v^ that enter 
the definition of A^. The obtained values are v\$ = 0.9019, Uu = 1-5010, z/ 23 = 1.0483, 
z/ 24 = 1.7447. 

After these considerations the free parameters are reduced to the five coupling constants 
A13, A 2 3, A14, A 2 4 and An. First of all we solved the imaginary- axis Eliashberg Equa- 
tions ([1]) and (J2J) (actually we continued them analytically on the real-axis by using the 
approximants Pade technique) and we fixed the free parameters in order to reproduce 
the gap values at low temperature. 

The large number of free parameters (five) may suggest that it is possible to 
find different sets that produce the same results. It is not so, as a matter of fact 
the predominantly interband character of the model drastically reduce the number of 
possible choices. 

At the beginning, in order to have the fewest number of free parameters, we set On 
to be the same of the antiferromagnetic SF (even if the intraband coupling can not be 
mediated by SF). In this case (see third row of table EJ the value of An is very large. 
We add only An because only the band 1 has a low Fermi energy and only in this band 
the Migdal's theorem breaks down. The effect of the vertex corrections [281 122] can be 
simulated by an effective coupling that is bigger than real coupling [301 El] • Certainly 
in the bands 2, 3 and 4, the phonon couplings are very small and therefore we have 
not considered the possibility to have A22, A33, A44 7^ and also, in a previous paper, 
we have demonstrated that, in these systems, the effect of a small phononic intraband 
coupling is negligible [24] . 

Then we considered for O n the typical phonon energies [7]. In this case (as 
reported in fourth and fifth row of table [2]) the value of An is 0.86-0.9, while 
the antiferromagnetic spin fluctuations contribution correspond to a moderate strong 
coupling regime {X s tot ~ 1.5). 

In fact, we have solved the Eliashberg equations in two other cases: in the first case 
we used as af-^Ffl) the calculated phonon density of states Gift) and the second case 
we considered the calculated total electro n-phonon spectral function a 2 F tot {fl) both 
appropriately scaled. The proper choice is the second one, but this spectral function 
should be in principle different for each band; since the first band shows peculiar 
characteristics the evaluation of the electron phonon coupling could be not so reliable. 
Then we decided to use the phonon spectra and let the coupling to be a free parameter. 
These spectral functions are shown in Figure [2l 

At this point there are no more free parameters. The critical temperature can be 
evaluated and it turns out to be very close to experimental one, T^ alc = 18.6 — 20.1 K. 
By solving real-axis Eliashberg equations we obtained the temperature dependence of 
the gaps (see Figure |3j) for the parameter sets reported in the third and fourth row of 
Table [2 

The multiband Eliashberg model developed above can also be used to explain the 
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Table 2. The first row shows the experimental data. The second row concerns the 
pure interband case (Xa = 0) while the last three include an intraband term (An ^ 0): 
A very large value appears in the first case (the third row) , a smaller one if the phonon 
spectral function G(fi) (fourth row) or the electron-phonon spectral function a 2 F(il) 
(fifth row) are considered. The critical temperatures are given in K and the gap values 
in meV. 
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Figure 2. The calculated phononic density of states G(f2) (red dashed line) and the 
calculated total electron-phonon spectral function a 2 F tot (£l) (black solid line). 



experimental data of temperature dependence of the upper critical magnetic field [19|. 
For the sake of completeness, we give here the linearized gap equations in the presence 
of magnetic field, for a superconductor in the clean limit. In the following, vpj is the 
Fermi velocity of the j — th band, and H c2 is the upper critical field: 
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Figure 3. Temperature dependence of the absolute gap values (lines) and 
experimental data (symbols) at 8K. The dark cyan solid (dashed) line represents the 
first gap, the orange solid (dashed) line the second one, the violet solid (dashed) line 
the third and the red solid (dashed) line the fourth, calculated with the parameters of 
fourth (third) row of Table EI 
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Figure 4. Experimental temperature dependence of the upper critical held (symbols) , 
and the relevant fitting curves (lines) obtained by solving the Eliashberg equations in 
magnetic field. Red circles and solid red (dashed dark blue) line for Hu c , black square 
symbols and solid black (dashed orange) line for H^ ab calculated with the parameters 
of fourth (third) row of Table [2j The dotted olive (H\\ c ) and navy (H\\ ab ) and the 
dashed-dotted magenta (H\\ c ) and pink (-ff|| a fc) lines are, respectively, the fourth and 
first case of Table [2] but where vfi isn't a free parameter but it is taken by Table [1] 



Here (3j = 7rif c2 f^ J /(2$o) and $o is the unit of magnetic flux. In these equations the 
four bare Fermi velocities vpj [32J are the input parameters. As the first band shows 
peculiar characteristics even in the calculation of the Fermi velocity can be present some 
corrections, then we decided to let the first Fermi velocities to be free parameters and 
we choose them to find the best fit of the experimental data [19] while the other values 
have been fixed to the values reported in Table [U Then vfi, in each case, is the only 



free parameter. The obtained values are: 



J F\ 



2.28 • 10 5 m/s and vf{ = 1.74 • 10 
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m/s, in the phonon case and Vp X = 2.79 • 10 5 m/s and vp^ = 2.14 • 10 5 m/s, if the spin 
fluctuation spectral function is considered. Figure H] shows the experimental data and 
the best theoretical curves (solid and dashed lines) obtained by solving the Eliashberg 
equations within the model discussed above. As can be seen, the results obtained in 
the two considered cases are almost indistinguishable and in very good agreement with 
the experimental data. In the Figure H] we shows also the curves calculated with vfx 
taken by Tabled] when An ^ (fourth case in Table |2j dotted line olive and navy) and 
when An — (first case in Table [21 dashed-dotted line magenta and pink). In both 
situations there is no agreement with the experimental data. The curves calculated in 
absence of the term An don't agree with the experimental data so we deduce that the 
higher value of vp\ ' c isn't produced by the presence of an intraband term (An ^ 0) but, 
probably, by the peculiar characteristics of band 1. However one must consider the fact 
that the Eliashberg equations are derived by assuming Migdal's theorem. In presence 
of an anomalous band dispersion as for band 1, the theory may partially break down. 
Allowing Vp^ c as a free parameter implicitly implies that we are "phenomeno logically" 
going beyond the first order contributions (i.e. now we cannot neglect the effects of the 
vertex corrections in the band 1). The break down of the Migdal's theorem leads to use 
effective values of An and Vpi different from real value because the framework of the 
theory is partially inadequate. 

To summarize we have constructed a phenomenological model of superconductivity 
for LiFeAs able to describe its critical temperature, the multigap structure measured by 
Umezawa and coworkers in Ref. [18] and other experimental observations. However this 
process was not straightforward, to be able to conjugate a spin fluctuation dominated 
pairing with the experimental gap structure we have been forced to introduce an 
intraband coupling that acts only on the first band. The presence of a phononic, purely 
intraband term seems to indicate an intrinsic incompatibility between the structure of 
the superconducting gaps as measured by Umezawa and coworkers in Ref. [18] and a 
purely spin-fluctuation mediated pairing. A possible explanation may be linked to the 
very low Fermi energy of the band for which vertex corrections [281 129] to the usual 
Migdal-Eliashberg theory may be relevant, and are expected to increase the stenght of 
the phononic pairing [301 E] • 

In conclusion, our calculations show that LiFeAs presents peculiar features with respect 
to other iron compounds and it cannot be explained within the framework of a pure 
interband spinfluctuation mediated superconductivity. 
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